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ABSTRACT 

Aims. The synthesis of the cosmic X-ray background (CXB) requires a large population of Compton-thick active 
galactic nuclei that have not been detected so far. We probe whether bright infrared galaxies could harbor a population 
of Compton-thick nuclei and if they could contribute significantly. 

Methods. We analyzed 112 Msec of INTEGRAL observations obtained on 613 galaxies from the IRAS Revised Bright 
Galaxy Sample. We derived the average hard X-ray (18-80 keV) emission of Seyfert and various non Seyfert galaxy 
subsamples to estimate their relative contribution to the locally emitted CXB. 

Results. The Seyfert 1 & 2 are detected at hard X-rays. None of the other galaxy subsamples were detected. ULIRGs are 
at least 5 times under-luminous at hard X-rays when compared to Seyferts. The upper limit obtained for the average 
non Seyfert galaxies is as low as 7 x 10"^'^ erg s~^ cm~'^. On average, these galaxies do not contain active nuclei brighter 
than lO*^ erg/s at hard X-rays. The total hard X-ray flux detected from the sample is 4.9 x 10~^ erg s~^ cm~^ (about 
1% of the CXB), and 64% of this originates in absorbed active nuclei. Local non-Seyfert galaxies contribute for less 
than 7% and do not harbor the Compton-thick nuclei assumed to synthesize the locally emitted CXB. 

Key words. Gamma-Rays: observations - X-rays: galaxies - X-rays: diffuse background ~ Infrared: galaxies - Galaxies: 
Seyfert 



1. Introduction 

While accreting super-massive bla ck-holes are believe d to 
power active galactic nuclei (AGN; Lynden-Bell|1969 1 , the 
apparent spectral energy distribution of these objects de- 
pends on many physical parameters, such as the black-hole 
mass, the accretion rate, the nature of the accretion flow, 
the inclination on the line of sight, and the intrinsic obscu- 
ration. 

The hard X-ray background, with a prominent maxi- 
mum around 30 keV, re quires that most AGNs are obscured 
( Setti fc Woltjer|[l989 1. Quantitatively, population synthe- 
sis models indicate that 60 to 80% of all AGNs should be 
obscur ed and that about 50% of these should be Compton- 
thick (Gilli et al. 2007), i.e. featuring absorbing column 
densities A^h > ~ 1-6 x 10^'* cm~^. Nuclear starbursts 
in the central few parsecs of a galaxy could p rovide the 
necessary obscuring material (Ballantyne '2008) to explain 
a large population of Compton-thick objects. These nuclear 
starbu rsts also link stellar formation and the black-h ole ac- 
tivity ( [Watabe et al]|2008| [Kawakatu fc Wada|[2008 1 . 

Various techniques have been used to study the rela- 
tion between starbursts and AGN activity, based on the 
study of th e infrared spectral energy distribution (Genzel 



et al.pMl or on c orrelations between infrared and X-ra y 
or radio ■detections ( [Barger et al.|2005[|Polletta et al.|2006| . 
Because hard X-ra y emission is an ubi quitous signature of 



AGN activity (Mushotzky et al. 1993 1 and as it is insen- 



sitive to photoelectric absorption, the correlation between 



infrared and hard X-ray fluxes is an interesting probe for 
the existence of highly absorbed AGNs. 

Even if a fraction of the AGN population features var- 
ious signatures of absorption - like an ab sence of broad 
emission lines, scattered emission (Schmidt et al. 2007) 



reddening ( [Walter fc Fink||1993[ [Wilkes et al.||2005|), or 



soft X-ray ]Dhotoelectric absorption (Tozzi et al. 2006 ) 
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recen t low sensitivity ha rd X-ray surveys (Paltani et al. 
2008| [Tueller et ar]|2008| ) have failed to detect the popula- 
tion of Compton-thick objects detected at low X-ray fluxes 
dRisaliti et al]|1999[ [Guainazzi et aL]|2005[ ). 

To improve the sensitivity of the current hard X-ray 
surveys by one order of magnitude, we stacked the hard 
X-ray signal measured by INTEGRAL for various types of 
infrared bright galaxies to check their contribution to the 
hard X-ray background. The sample of sources is presented 
in section |2l the INTEGRAL data and analysis is described 
in section l3| and their correlation with infrared data is dis- 
cussed in section m 



2. Source sample 

In 1983, the Infra-Red Astronomical Satellite (IRAS) sur- 
veyed approximately 96% of the sky in four broad wave- 
length bands centered on 12, 25, 60 and 100 /.tm. The IRAS 
Point Source Catalog (PSC) and the IRAS Faint Source 
Catalog (FSC) were produced hsting 245889 and 173044 in- 
frared sources respectively (the PSC extends to lower galac- 
tic latitude than the FSC). IRAS discovered the luminous 
and ultra-luminous infrared galaxies (LIRGs and ULIRGs) 
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Table 1. Sources from the RBGS excluded from the anal- 
ysis and reason for their exclusion. 



Sources 


Reason 


NGC 2339, IRAS 05405+0035, 
NGC 2342, IRAS 05442+1732, 
IRAS 05223+1908 


Close to Crab 


NGC 2993, NGC 3583, NGC 4402, 
NGC 4303, NGC 6215, NGC 7590, 
NGC 7599, IRAS 18293-3413 


Close to other bright 
sources 


LMC, SMC 


Extended 


NGC 6240 


Detected, complex 



in the local Universe with infrared luminosities larger than 
i Q re spectively. 
pOBj ) 



IQii L^. and lO^^ L, 



Sanders et al. 



used the IRAS FSC and PSC 
to constitute the IRAS Revised Bright Galaxy Sample 
(RBGS). The RBGS is a complete flux-limited sample of 
all extragalactic objects brighter than 5.24 Jy at 60 /im, 
covering the entire sky surveyed by IRAS at galactic lati- 
tudes |b| > 5°. The median redshift of the RBGS galaxies is 
z=0.0082 and their bolometric luminosity function breaks 
at ~ lO^^-^L© from a slope of -0.6 to -2.2. 

Among the 629 objects of the RBGS, we excluded 
(Table [T|) the SMC and LMC because they are extended 
for the imager on board INTEGRAL, 13 objects because 
of their proximity with other sources detected at hard X- 
rays and NGC 6240, one of the rare object hosting two 
massive black-holes. NGC 6240 is an ULIRG and includes 
a starburst, a Com pton-thick and a normal Seyfert nuclei 

I Risaliti et al.|2006 1 . The origin of the hard X-ray emission 
Vignati et al.||1999| ) detected by INTEGRAL at a level of 
mCrab is therefore unclear. 
The resulting sample of 613 objects was classified in 
six sub-samples according to their object type defined by 
SIMBAeQ 93 Seyfert galaxies (including 14 Seyfert 1, 64 
Seyfert 2 and 15 Seyfert of unclassified or complex types) 
and 520 "non Seyfert" galaxies, missing clear signature 
for an active nucleus, including 12 ULIRGs (Leo/^m/L© > 
10^2), 144 LIRGs (10^2 > Leo^nn/L© > 10") and 41 Liners. 
Among the 14 Seyfert 1 galaxies of our sample, several ob- 
jects have in fact intermediate Seyfert types. NGC 1386, 
NGC 2992, NGC 5033 are Seyfert 1.9 bu t their X-ray 



spectra do not show sig nificant absorption (Dadina 



Terashima et al. 


19991 


Seyfert 1.5 with complex 



2008 



NGC 3227 and NGC 4l5l~are 
md variable absorption ('Gondoin' 
et al.|2003|[Zdziarski et al.,2002^ . The division of the Seyfert 
sample according to their Seyfert type is only an approx- 
imation of their ab sorption properties in the X-rays (e.g. 
Bassani et aT]|2006| ). 



3. INTEGRAL data and analysis 

All pubHc data obta ined by the INTEGRAL soft 7-ray 



image r IBIS/ISGRI (Ubertini et al. 2003 



2003 1 until January 2007 are considered in this study. The 



Lebrun et al. 



pointings, with typical exposures of (1 — 3) x 
ing at least one source of the sample in its field of view and 
with an effective exposure longer than 120 s were selected, 
spanning times between Dec. 30, 2002 (revolution 26) and 
Jan. 17, 2007 (revolution 520). 

The ISGRI data were reduced using the INTEGRAL 
Offline Scientific Analysis softwar^ version 7.0 pub- 
l icly released by the IN TEGRAL Science Data Centre 
(Courvoisier et al. 2003). Individual sky images for each 



INTEGRAL imager has a large field of view of 29° square 
with a spatial resolution of 12 arcmin. Several thousands 



pointing have been produced in a broad energy band of 18- 
80 keV. All known ISGRI sources, with significance above 
5(7, have been used for image cleaning. The position of these 
sources was fixed to catalogued values when known accu- 
rately. 

To measure the average hard X-ray properties of the 
selected galaxy subsamples we have built mosaic images 
of all subsample sources modifying the coordinate system 
of each individual image such that the coordinates of each 
subsample sources were set to an arbitrary fixed position 
{a—0, 5=0). The resulting image provides a stack of all 
ISGRI data, available on all sources of interest. 

To minimize the systematics in the mosaic images, 
we have excluded 1455 individual sky images (many ob- 
tained before revolution 38, when the IBIS bottom anti- 
coincidence was reconfigured) which had background fiuc- 
tuation rms larger than 1.1 cr in the significance image and 
779 images with minimum significance lower than — 5.5(7. 
36030 images were finally included in the processing (some 
of them many times, when including several sources). 

The 500x500-pixels mosaic images were built in equato- 
rial coordinates centerred on (a=0, (5=0) with a tangential 
projection using a factor two over-sampling when compared 
to the individual input sky images; this results in a pixel 
size of 2.4 arcmin in the center of the mosaic. The photo- 
metric integrity and accurate astrometry are obtained by 
calculating the intersection between input and output pix- 
els and weighting count rates according to the overlapping 
area. A composite of the central part of three mosaic images 
is featured in Fig. [T] 

The effective exposure obtained at the center of the mo- 
saics are between 1.7 and 96.0 Ms for the various galaxy 
subsamples. To investigate the quality of the mosaics, we 
studied the distribution of the pixels significance, which is 
expected to be Gaussian with average and dispersion 1. 
We fitted a Gaussian distribution to the part of the his- 
togram with (7 ^ 3 to avoid the strong positive tail due 
to real sources. While the width of the pixel distributions 
(« 1.2) are not unusual the centroid of the distribution is 
significantly displaced towards positive values for the sam- 
ple of "non Seyfert" galaxies where the effective exposure 
reaches almost 100 Ms. Our analysis shows that the indi- 
vidual ISGRI images appear slightly positive on average 
(0.004 ct/s) which could not be detected with typical ex- 
posures up to several Ms. Even if the exact cause of this 
effect is not known (background subtraction ?), this effect 
can be corrected for by subtracting a constant value to the 
mosaic image. 

A positive average signal is clearly detected in the mo- 
saics built for the three Seyfert subsamples. No average 
excess is detected for the other galaxy subsamples. As ex- 
pected, the position {a=0, S=0) is well within the 90% 
error box corresponding to the excesses found for the three 
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Fig. 1. 18-80 keV stacked mosaic significance images for 
three of the various galaxy subsamples. From top to bot- 
tom: Seyfert 1, Seyfert 2 and "non Seyfert" subsamples. 
The color map spans significances between 1 and 6. The 
color map of the bottom image has been slightly modified 
to account for the systematics, re-normalizing the centroid 
of the background distribution to a value of 0. 



Seyfert subsamples and the la width of these excesses are 
between 5.7 and 6.0 arcmin, in good agreement with the 
ISGRI point spread function. 

Table [2] lists the average signal (or 5a upper limits) 
extracted at hard X-rays from each ISGRI mosaic with 
mosaic_spec (corrected for the higher sky background sys- 
tematics obtained for the large "non Seyfert" galaxy sub- 
sample) to gether with the average flux measured by IRAS 
at 60 /Ltm (Sanders et al. 2003) for each subsample. The 
hard X-ray fluxes were derived from the count rates using 
a spectral index of F = 1.8. The hard X-ray detections ob- 
tained for the three Seyfert subsamples are more than one 
order of magnitude brighter than the upper-limits obtained 
for the other galaxy subsamples with comparable effective 
exposures. The hard X-ray upper limit obtained for the 



Table 2. Average signals (or 5a upper limits) detected for 
the various subsamples at hard X-rays (18-80 keV) and 
at 60 fim. Effective hard X-ray exposure, count rates and 
integrated flux are provided together with the average 60 
/Ltm flux. The hard X-ray detection signiflcances for the first 
four samples are 119, 23, 122 and 147 respectively. 



Sample 


Srcs 


Exp. 


Ct rate 


Fl8-80kcV 








Ms 


ct/s 


10^^" Cgs 


Jy 


Seyfert 1 


14 


1.7 


2.25±0.02 


1.28 


18.3 


SI - 4151^ 


13 


1.4 


0.51±0.02 


0.29 


19.3 


Seyfert 2 


64 


12.1 


0.87±0.01 


0.49 


31.1 


Seyfert 


93 


16.1 


1.04±0.01 


0.59 


28.9 


ULIRGs 


12 


4.1 


< 0.06 


< 0.034 


10.2 


Liners 


41 


7.5 


< 0.045 


< 0.026 


29.6 


LIRGs 


144 


26.7 


< 0.02 


< 0.011 


12.9 


non Seyfert 


520 


96.0 


< 0.013 


< 0.0074 


19.1 



^includes all Seyfert 1 galaxies excepting NGC 4151. 



"non Seyfert" subsample is a hundred times less than the 
flux detected in average from the Seyfert galaxies. 

At hard X-rays, NGC 4151 dominates the Seyfert 1 sub- 
sample. Excluding it reduces the average count rate from 
2.25 to 0.51. The average count rate obtained for the small 
Seyfert 1 subsample is signiflcantly affected by NGC 4151, 
contrasting with the average 60 /im flux that is almost not 
affected. 



4. Discussion 

4.1. Average spectral energy distribution 

At 60 /im most of the the emission is believed to originate 
from molecular clouds or diffuse dust reprocessing stellar 
light or from nuclear activity for luminous AGNs. Figure 
|2] presents the re-normalized average spectral energy dis- 
tributions of each subsamples, together with the average 
i nfrared-optical s pectra obtained for a sample of 30 ULIRGs 
( Vega et al.|2008 ) for illustration. The spectra of the galaxy 
subsamples were re-normalized at 60 /im to the average ob- 
tained for the 30 ULIRGs. 

The ratios between the flux densities at hard X-rays and 
at 60 /im are 12% for Seyfert 1, 3% for Seyfert 2 (and for 
Seyfert 1 excluding NGC 4151), below 0.6% for ULIRGs, 
below 0.1% for the LIRGs and Liners, and flnally below 
0.07% for the "non Seyfert" galaxies all together. 

From the hard X-ray point of view we flnd a clear sep- 
aration between the average Seyfert dominated and non 
Seyfert galaxies, the latter being not detected. As the hard 
X-rays are not much affected by moderate absorption this 
difference could be intrinsic to the source samples and/or 
signature of Compton-thick absorption in the galaxy sub- 
samples lacking Seyfert signatures. 

The ratio between the flux densities obtained at hard X- 
rays from high-mass X -ray binaries and in the f ar infrared 
can be estimated from Grimm et al. ( 2003 1 and Kennicutt 
( 1998 ) and is not larger than 0.03%. This is two orders of 



magnitude lower than the hard X-ray signal obtained for 
the Seyfert subsamples and still below the upper-limits ob- 
tained for the "non Seyfert" galaxies. Our upper-limits are 
therefore consistant with the hard X-ray emission expected 
from stellar activity. The ratio of the 60 /im to 18-80 keV 
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Fig. 2. Average spectral energy distribution (SED) for the 
various galaxy subsamples, renormalized to an arbitrarily 
fixed 60 flux. The average ULIRGs infrared-optical SED 
is represented as a thin line for illustration. The hard X- 
ray spectral shape assumes a spectral slope of F = 1.8. The 
long dash curve shows the relative hard X-ray contribution 
of Seyfert 1 galaxies, excluding NGC 4151. 



flux densities in NGC 1068 {k. 500; [Matt et al.||1997[ ) is 
also smaller than the one obtained on average for the "non 
Seyfert" galaxies, in agreement with a lack of reflected hard 
X-rays in the latter. 

The Seyfert 2 are 4.5 times more numerous than the 
Seyfert 1 in the RGB S in agreement with the r esults of |La| 



Franca et al. (20051 and Akylas et al. (2006) for a sam- 



pie dominated by low luminosity sources. The hard X-ray 
emission of moderately absorbed objects is, at maximum, 
a factor of two lo wer than that of unabsorbed objects (e.g. 
Gilli et al. 2007). When compared to their 60 [im emis- 
sion, local Seyfert 1 appears on average 5 times brighter at 
hard X-rays than Seyfert 2. As most of this difference can 
be explained by the presence of NGC 4151, the sample is 
too small to elaborate on a possible genuine difference be- 
tween the average unabsorbed spectral energy distributions 
of Seyfert 1 & 2. 

The hard X-ray upper-limit obtained for the ULIRGs 
indicates that, on average, these sources emit less than 20% 
of the average emission detected in Seyferts when normal- 
ized to their infrared emission. This is consistent with the 
fact that less than 30% of the ULIRGs harbor a significant 
(> 10%) AGN contribu tion in the far infrared ( Armus et al. 
2007 Vega et al.|2008 1 and with the pre dominance of their 
starb urst emission even at soft X-rays (Franceschini et al. 
2003 1. Note that as the number of nearby ULIRGs is small. 



our upper-limits are an order of magnitude above the ex- 
trapolation at hard X-rays of the XMM- Newton detections. 

Active nuclei absorbed by Compton-thick material in 
ULIRGs or LIRGs could, in principle, remain completely 
hidden at hard X-rays. To accomodate our hard X-ray 
upper-limits and assuming that the active nuclei have the 
same average intrinsic spectral energy distribution than ob- 



tained for our sample of Seyfert galaxies, no more than 3 
and 20% of the hard X-ray emission should be transmit- 
ted on average in LIRGs and ULIRGs respectively. For the 
ULIRGs, this means that active nuclei, if they exist, should 
be obscured by Nh > lO^''^ cm~^ and that no more than a 
few of them in the 12 sources sample should remain unab- 
sorbed. For the LIRGs, the active nuclei, if present, should 
be obscured by Nh > 5 x 10^^ cm~^ and no more than a 
few out of 144 should be unobscured. 

The "non Seyfert" galaxies are, on average, very faint 
at hard X-rays. The upper-limit we have obtained indicates 
that their hard X-ray emission is not driven by active nuclei 
brighter than Lx ^ 10**^ erg/s on average. The presence of 
bright active nuclei absorbed by Nh > 5 x 10^** cm~^ with 
no more than 15 unobscured ones out of 520 sources can 
however not be ruled out. 

The hard X-ray upper-limit obtained for the Liner 
galaxies is in agreement with active nuclei or starbursts 
of low X-ray luminosity 10^° erg /s) reported in these 
objects ( Gonzalez-Martin et al.|[2006 l. 



4.2. Contributions to the hard X-ray background 

We can estimate the relative contributions of our subsam- 
ples of galaxies to the locally emitted hard X-ray back- 
ground (CXB). The cumulative hard X-ray flux of the 
Seyfert galaxies is 1.8 x 10~^ and 3.1 x 10^^ erg s~^ cm~^ 
for the type 1 a nd 2 respectively (shghtly less than 1% of 
the total CXB; [Gilli et al 2007 and references therein), 
e hard . 



The cumulative 



X-ray flux of the non Seyfert galaxies 

1 „v„-2 



is < 0.38 X 10"" erg s~" cm 

Assuming that Seyfert 1 and 2 are optically thin and, 
respectively, mildly absorbed (which is not completely ac- 
curate in gen eral; e.g. |de Rosa et al.|[2008{ [Brightman fc 
Nandra 2008), the fraction of the locally emitted CXB 



emitted by absorbed objects (Nh > 10^^ cm^^) in our 
sample, is 64% and the fraction emitted by potentially 
Compton-thick non Seyfert objects is < 7% (5(t). 

These fractions can be compared to the predictions of 



et al. 
tool 



CXB synthe sis models and in particular to these of Gilli 
(2007), which could be obtained through the on-line 
'OMR^ Our sample is flux limited in the infrared, 
the corresponding X-ray flux limit for an individual Seyfert 
galaxy is on the order of 10^^^ erg s~^ cm~^. For such a 
flux limited sample the CXB synthesis model indicates that 
absorbed and Compton-thick sources should contribute for 
65% and, respectively, 18% of the locally emitted CXB at 
30 keV. 

The fraction of the locally emitted CXB to be accounted 
for by absorbed sources agrees with the results of our obser- 
vations. The fraction expected to be emitted by Compton- 
thick sources is however signiflcantly larger than the upper- 
limits obtained for non-Seyfert galaxies, which are therefore 
unlikely to harbor Compton-thick active nuclei contribut- 
ing significantly to the locally emitted CXB. 

5. Conclusions 

We have derived average hard X-ray fluxes for sub- 
samples of Seyfert and non Seyfert galaxies stacking 
INTEGRAL/ISGRI images of 613 members of the IRAS 
Revised Bright Galaxy Sample. 



http: / /www. bo. astro. it /~gilli /counts. html 



R. Walter and N. Cabral: The hard X-ray view of bright infrared galaxies 



5 



The Scyfcrt 1 & 2 subsarnplcs arc clearly detected at 
hard X-rays while none of the other galaxy subsamples is 
detected. When compared to their 60 /xm emission, local 
Seyfert 1 are on average 5 times brighter at hard X-rays 
than Seyfert 2. This can be explained by the presence of 
NGC 4151 dominating the small subsample of Seyfert 1. 

The upper-limit obtained on the ULIRGs indicates that 
on average these objects emit less than 20% of the average 
Scyfcrt emission at hard X-rays. This is consistant with pre- 
vious results suggesting that the majority of these objects 
do not harbor significant AGN activity. 

On average, "non Seyfert" galaxies do not contain any 
active nuclei brighter than 10*^ erg/s at hard X-rays. If 
they harbor Compton-thick AGNs, these should be highly 
absorbed (> 5 x 10^"^ cm~^) and no more than a few % of 
them could remain unabsorbed. 

The RGBS accounts for less than 1% of the CXB. 64% 
of this is produced by absorbed (> lO'^^ cm^^) sources, 
in very good agreement with CXB synthesis models. Non 
Seyfert galaxies contribute for less than 7% {5a) and have 
not been detected statistically. These galaxies are unlikely 
to harbor the Compton-thick AGNs required to synthesize 
the locally emitted CXB. Local Compton-thick AGNs, if 
they exist as a population, should therefore be looked for 
among Seyfert 2 or infrared weak galaxies. 
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